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Abstract 

In this paper we answer the question of when circulant quantum spin networks with nearest-neighbor 
couphngs can give perfect state transfer. The network is described by a circulant graph G, which is 
characterized by its circulant adjacency matrix A. Formally, we say that there exists a perfect state 
transfer (PST) between vertices a,b € V{G) if |-F'(T)ab| — 1, for some positive real number r, where 
F{t) — exp{iAt). Saxena, Severini and Shparlinski (International Journal of Quantum Information 5 
(2007), 417-430) proved that \F{T)aa\ = 1 for some a G V{G) and r e M+ if and only if all eigenvalues 
of G are integer (that is, the graph is integral). The integral circulant graph ICG„(-D) has the vertex 
set Zn — {0, 1, 2, . . . , n — 1} and vertices a and b are adjacent if gcd(a — b,n) G D, where D C {d : d \ 
n, 1 < d < n}. These graphs are highly symmetric and have important applications in chemical graph 
theory.^e show that ICG„{D) has PST if and only if n e 4N and D = U D2 U 2D2 U 4D2 U {n/2''}, 
where D3 = {d & D \ n/d & 8N}, D2 = {d £ D \ n/d & 8N + 4} \ {n/4} and a G {1, 2}. We have thus 
answered the question of complete characterization of perfect state transfer in integral circulant graphs 
raised in Quantum Information and Computation, Vol. 10, No. 3&4 (2010) 0325-0342 by Angeles-Canul 
et al. Furthermore, we also calculate perfect quantum communication distance (distance between vertices 
where PST occurs) and describe the spectra of integral circulant graphs having PST. We conclude by 
giving a closed form expression calculating the number of integral circulant graphs of a given order having 
PST. 
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1 Introduction 



The transfer of a quantum state from one location to another is a crucial ingredient for many quantum 
information processing protocols. There are various physical systems that can serve as quantum channels, 
one of them being a quantum spin network. These networks consist of n qubits where some pairs of qubits 
are coupled via XY-interaction. The perfect transfer of quantum states from one qubit to another in such 
networks was first considered in [7]. There are two special qubits A and B representing the input and output 
qubit, respectively. The transfer is implemented by setting the qubit A in a prescribed quantum state and 
by retrieving the state from the output qubit B after some time. The transfer is called perfect state transfer 
(transfer with unit fidelity) if the initial state of the qubit A and the final state of the qubit B are equal up 
to a local phase rotation. 

Every quantum spin network with fixed nearest-neighbor couplings is uniquely described by an undirected 
graph G on a vertex set V{G) ={1,2,..., n}. The edges of the graph G specify which qubits are coupled. 
In other words, there is an edge between vertices i and j if i-th and j-th qubit are coupled. 

In [7] a simple XY coupling is considered such that the Hamiltonian of the system has the form 
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and crf,cr^ and cr? are Pauli matrices acting on i-th qubit. The standard basis chosen for an individual 
qubit is {|0), |1)} and it is assumed that all spins initially point down (|0)) along the prescribed z axis. In 
other words, the initial state of the network is |0) = |OyiO . . . 00^). This is an eigenstate of Hamiltonian Hq 
corresponding to zero energy. The Hilbert space Tic associated with a network is spanned by the vectors 
\e1e2 ■ ■ ■ en) where E {0, 1} and, therefore, its dimension is 2". 

The process of transmitting a quantum state from A to B begins with the creation of the initial state 
a|OyiO . . . 00b) + /3|1aO . . . 00^) of the network. Since |0) is a zero-energy eigenstate of Hq, the coefficient a 
will not change in time. Since the operator of total z component of the spin a^^f- = Y17=i commutes with 
Hg, state 1 1^0 . . . 00^) must evolve into a superposition of the states |i) = |0 . . . OljO, . . . , 0) for i = 1, . . . , n. 
Denote by Sq the subspace of Hg spanned by the vectors i = 1, . . . , n. Hence, the initial state of network 
evolves in time t into the state 

n 

a\0) + Y,Pi{t)\i) eSG- 

i=l 

The previous equation shows that system dynamics is completely determined by the evolution in n-dimensional 
space Sg- The restriction of the Hamiltonian Hg to the subspace Sg is an n x n matrix identical to the 
adjacency matrix Ag of the graph G. 

Thus, the time evolution operator can be written in the form F{t) = exp(MGt). The matrix exponential 
exp(M) is defined as usual 

+00 ^ 

exp(M) = -M^ 

n=0 

Perfect state transfer (PST) between different vertices (qubits) a and 6 (1 < a, 6 < n) is obtained in time 
r, if {a\F{t)\b) = \F[T)ab\ = 1- The graph (network) is periodic at a if \F[T)aa\ = 1 for some r. A graph is 
periodic if it is periodic at each vertex a. 

The existence of PST for some network topologies has already been considered in the literature. For 
example, Christandl et al. [8] proved that PST occurs in paths of length one and two between their 
end-vertices and also in Cartesian powers of these graphs between vertices at maximal distance. In the 
recent paper Godsil constructed a class of distance-regular graphs of diameter three, with PST. Some 
properties of quantum dynamics on circulant graphs were studied in [J. Saxena, Severini and Shparlinski 
|18j considered circulant graphs as potential candidates for modeling quantum spin networks having PST. 
They show that a circulant graph is periodic if and only if all eigenvalues of the graph are integers (i.e. 
graph is integral). Since periodicity is a necessary condition for PST existence |18] . circulant graphs having 
PST must be integral circulant graphs. A simple and general characterization of the existence of PST in 
an integral circulant graph, in terms of its eigenvalues, was given by Basic, Petkovic and Stevanovic in [3]. 
Furthermore, it was shown that for odd number of vertices, there is no PST, and that among the class of 
unitary Cayley graphs a subclass of integral circulant graphs, only K2 (the complete graph with two nodes) 
and C4 (the cycle of length four) have PST. In the recent paper [T7], it was proven that there exists an 
integral circulant graph with n vertices having PST if and only if 4 | n. Several classes of integral circulant 
graphs having PST were found as well and several others in jl5j. 

In all known classes of graphs having PST perfect quantum communication distances (i.e. the distances 
between vertices where PST occurs) are considerably small compared to the order of the graph. One idea for 
the distance enlargement, is to consider networks with fixed but different couplings between qubits. These 
networks correspond to graphs with weighted adjacency matrices. For example, in [TJIH] the authors showed 
that PST can be achieved over arbitrarily long distances in a weighted linear paths. Many recent papers 
have proposed such an approach [HI [151 ES] ■ 

Studying PST in integral circulant graphs can also be interpreted as a contribution to the spectral theory 
of integral graphs. These graphs are highly symmetric and have some remarkable properties connecting 
graph theory and number theory. The term 'integral circulant graph' first appears in the work of So [20j . 
where a nice characterization of these graphs in terms of their symbol set is given. The upper bounds on the 
number of vertices and the diameter of integral circulant graphs were given in |18j . Furthermore, Stevanovic, 
Petkovic and Basic [19] improved the upper bound. Various other properties of unitary Cayley graphs were 
recently investigated. For example, Berrizbeitia and Giudici [6] and Fuchs [9] established the lower and 
upper bound on the size of the longest induced cycle. Klotz and Sander jl3] determined the diameter, clique 
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number, chromatic number and eigenvalues of unitary Cayley graphs. Basic and Ihc [2] calculated the clique 
number of integral circulant graphs with exactly one and two divisors and also provided an inequality for 
the general case. 

In this paper we proceed with the study of circulant networks supporting PST initiated in [3l[31 [T5l[T71 [T8]. 
First we give some properties of the spectra of ICG,i(Di) where Di = {d £ D \ 4^ n/d}. In that subsection 
we present some preliminary results which are used in the sequel of the section. We show that ICG„(Z?) 
has PST if and only if n e m and D = D^U D2U 2D2 U AD2 U {n/l"}, where D3 = {d e D \ n/d e 8N}, 
D2 = {d £ D \ n/d G 8N + 4} \ {n/4} and a G {1,2}. Using the previous result we prove that perfect 
quantum communication distance is either equal to one if n/2 G D or equal to two if n/4 G D. Moreover, 
we describe the spectra of integral circulant graphs having PST. The paper is concluded with a formula for 
the number of lCGn{D) having PST as a function of the number of vertices n. This results answer questions 
posed in [lailTKIB]. 



2 Integral circulant graphs 

The circulant graph G{n; S) is a graph on vertices Z„ = {0, 1, . . . , n — 1} such that each vertex i is adjacent 
to vertices i +n s for all s G S. The set S C Z„ is called the symbol of the graph G(n; S) and +„ denotes 
addition modulo n. Note that the degree of G{n; S) is #5. A graph is integral if all its eigenvalues are 
integers. Wasin So has characterized integral circulant graphs [20j in the following theorem: 

Theorem 1 [20] A circulant graph G{n; S) is integral if and only if 

S=\J Gn{d), 

for some set of divisors D C Z)„. Here = {k : gcd(/c,n) = d, 1 < k < n — 1}, and Dn is the set of 

all divisors of n, different from n. 

Therefore an integral circulant graph (ICG) G{n] S) is defined by its order n and the set of divisors 
D. Such graphs are also known as gcd-graphs (see for example [13]). An integral circulant graph with n 
vertices, defined by the set of divisors D C Z)„ will be denoted by ICG„(-D). Prom Theorem [T] we have that 
the degree of an integral circulant graph is degICG„(D) = XldeZ) ^{n/d). Here ^p{n) denotes the Euler-phi 
function [12] . 

The eigenvalues and eigenvectors of ICG„(D) are given in [18] as 

A,=^a;^^ ^, = [la;>2s...^(n-i).]^ (1) 

ses 

where ojn = exp(i27r/n) is the n-th root of unity. Denote by c{n,j) the following expression 

, . , , , ip(n) n 
c(.,n) = M^„,);^, ^r^. = ^^y (2) 

where ji is the Mobius function defined as 



{1, ifn = l 
0, if n is not square-free (3) 
(—1)^, if n is product of k distinct prime numbers. 

The expression c{j,n) is known as the Ramanujan function ([121 P- 55]). Eigenvalues Xj can be expressed 
in terms of the Ramanujan function as follows ([13j, Theorem 16) 

A, = ^c(i,n/d). (4) 

Let us observe that the Ramanujan function has the following basic properties which we will make use of in 
the paper. 



3 



Proposition 2 For any positive integers n,j and d such that d \ n, holds 

c(0,n) = ip{n), (5) 

c(l,n) = ^(n), (6) 

{li{n), nG2N + l 

/x(n/2), nG4N + 2 (7) 

2/x(n/2), n G 4N 

c{n 2,n d) = < ) , r -, (8) 

Proof. Directly using relation □ 

The integral circulant graph ICGn(-D) is connected if and only if gcd(n, di, . . . , d^) = 1 where D = 
{di, . . . , dk}- In the rest of the paper we will only consider connected integral circulant graphs. 

3 Perfect state transfer 

Let G be an undirected graph and denote by Aq its adjacency matrix. Let F{t) = exp{iAGt). There is a 
perfect state transfer (PST) in graph G [TJ \TT\ [T8] if there are distinct vertices a and b and a positive real 
number t such that \F{t)ab\ = 1- 

Let Ao, A2, . . . , A„_i be the eigenvalues (not necessarily distinct) of Aq and uq, ui, . . . , Un-i be the cor- 
responding normalized eigenvectors. We use spectral decomposition of the real symmetric matrix Aq (see 
for example [lOj (Theorem 5.5.1) for more details). The matrix function F{t) can be represented as 

n-1 

F{t) = Y,^M-^>^kt)ukul. (9) 

k=0 

Now let G = ICG„(D) be an integral circulant graph. By simple calculation and using ([T]), we see that 
ll^^fcll = hence Uk = Vk/y/n. Expression ([9]) now becomes 



F{t) = -y^exp{iXkt)vkvl. 

k=0 

In particular, from the last expression and ([1]) it directly follows 

n-1 



F{t)ab = -y] exp(iAfct)a;^('^- 



n 

k=0 



This expression is given in [18j (Proposition 1). Finally, our goal is to check whether there exist distinct 
integers a, 6 G Z„ and a positive real number t such that \F{t)ab\ = 1, i-e. 



^ n— 1 

-^exp(iA,t)a;^(- 



k=0 



1. (10) 



Since the left-hand side of (|1U|) depends on a and b only as a function of a — 6 we can, without any loss 
of generality, assume that b = 0. Therefore, throughout the paper we consider the existence of PST only 
between vertices o and 0. 

We restate some results proved in [4J. These results establish necessary and sufficient conditions for (jlOp . 



Theorem 3 [4] There exists PST in ICG„(Z)) between vertices a and if and only if there are integers p 
and q such that gcd(p, q) = 1 and 

^(A,+i-A,) + - GZ, (11) 
q n 

for all j = 0, . . . ,n — 2. 
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Theorem 4 [4] There is no PST in ICGn(-D) if n/d is odd for every d £ D. For n even, if there exists 
PST in ICG„(-D) between vertices a and 0, then a = n/2. 

According to Theorem HI PST may exist in ICG„(L') only between vertices n/2 and (i.e., between h 
and n/2 + b as mentioned in [18]). Hence we will avoid referring to the input and output vertex and will 
just say that there exists PST in lCGn{D). 

The next corollary is derived from Theorem [3] and is further used as the criterion for a nonexistence of 
PST. 

Corollary 5 [4| // Aj = for some j = 0, . . . ,n — 2 then there is no PST in ICGn(D). 

For a given prime number p and an integer n G nJ^, denote by Sp{n) the maximal number a such that 
I n if n G N, and Sp{0) = +oo. The following result was proven in [4] and is further used as a criterion 
for the existence of PST. 

Lemma 6 [4] There exists PST in ICG„(-D), if and only if there exists a number m S No such that the 
following holds for all j = 0,1, ... ,n — 2 

S2{Xj+i - A,) = m. (12) 
The following corollary follows directly from Lemma [6l 

Corollary 7 Let ICG„(D) have PST. One of the following two statements must hold 

1. Xj = Xj+i (mod 2) for every < j < n — 1 (i.e., all eigenvalues Xj have the same parity). 

2. Xj = Aj+i + 1 (mod 2) for every < j < n — 1 (i.e., Xj are alternatively odd and even). 

We end this section with the following result concerning nonexistence of PST in ICG„(D), where n G 
4N + 2. 

Theorem 8 ^17j There is no PST in ICG„(D) for an arbitrary set of divisors D for n G 4N + 2. 

4 Integral circulant graphs having PST 

Let ICGn{D) be an arbitrary integral circulant graph. We define sets Di C D for < i < ^, where I = S2{n), 
in the following way 

Di = {d€ D\ S2{n/d) = i]. 

For simplicity of notation we also define sets Di, D3 C to be Di = Dq U Di and D3 = Ui>3 Dj. 
Let us introduce the notation kD for the set {kd \ d G D} for a positive integer k. 

4.1 Spectrum of the integral circulant graph ICG„(Di) 

In this section we deal with integral circulant graphs ICG„(-D) such that for each divisor d £ D it holds 
that 41^. In the rest of the section we will consider only such classes of graphs unless otherwise stated. 
In Lemmas [91 and fTOl we present some properties of the Ramanujan function. 

Throughout the section, we let n = 2"0p"^p2^ • . . . • , where pi < p2 < ■ ■ ■ < Pk are distinct primes, 
and ^ 1 for 1 < i < k and ao ^ 0. 

Lemma 9 For n > 2 it holds that c{j,n) G 2N + I if and only if A \ n and j = Pi^~^ • • ■p'^''~^J for some 
integer J such that gcd( J,n) G {1,2}. 

^No = N U {0} 
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Proof. 



(=^:) Suppose that c{j,n) is an odd integer. Since c{j,n) = fi{tn,j)^in-)/(p{tn,j), it holds that fJ-itnj) = ±1, 
i.e. tnj is square- free and ip{n) / ip{tn,j) is an odd integer. 

Suppose that for some odd pi it holds that pi f tnj- Let n' = n/p°'' . Since tnj \ n' and so ip{tnj) \ ^{n') 
we obtain that 

The last equation implies that c(j, n) is even since Pi — 1 is even. This is a contradiction and we can conclude 
that Pi I tn,j for every 1 < i < k. 

Now we have that ^{tnj) = {pi — I) ■ ■ ■ {pk — 1) and thus 

c(i,n) =2°»-ip°i-ip°2-i...p-.-i. 

Since c(j, n) is odd it holds that < ao < 1 or equivalently 4 f n. 

If n G 2N + 1 it must hold that tnj = Pi - ■ -Pk since tn,j is square-free. If n G 4N + 2, we have two 
possibilities for tn,j-, tnj = Pi - • - Pk or tnj = 2pi ■ ■ - pk depending on the parity of j. 

Furthermore, using n = gcd{n, j)tn,j we obtain that gcd(n, j) = • • (tnj and n have the same 

parity) or gcd(n,j) = 2p"^~^ • • -p^*"^ (otherwise). This implication of the lemma is now straightforward. 

(<;=:) Since gcd(n,j) = p"^^^^ ■ ■ ■ p'^''~^ gcd{J,n) and gcd(J, n) = {1,2}, it holds that tnj = Pi---pk or 
tn,j = 2pi ■ ■ -pk. In either case it holds that ip{tnj) = {pi — 1) ■ ■ ■ [p^ — 1). Now since 

c(j>) = ^^i^n,j)^^^ = ^vT~^ ■ • 
we conclude that c(j, n) G 2N + 1. □ 



Lemma 10 Let d he an arbitrary divisor of n such that n/d ^ 2N + 1 and < j < n — 1 he an arbitrary 
integer, then c{j,n/d) = —c{j,2n/d) for j G 2N + 1 and c{j,n/d) = c{j,2n/d) for j G 2N. 

Proof. As n/d G 2N + 1 we conclude that Lp{2n/d) = ip{n/d). 
Suppose that j G 2N + 1. Then gcd(2n/d, j) = gcd{n/d, j) and 

t '^^ 2 ^ 2t 

'^ri/d,j dgcd{2n/d,j) dgcd{n/d, j) "/'^'J' 

Furthermore, it holds that ^{t2n/d,j) = ^{in/dj) since ^n/dj is odd. Also we have that t2n/d,j is square-free 
if and only if tn/d,j is square-free, and ^J'{t2njd,j) = ~M(*n/d,j)- Now we can directly conclude that 

/ • r, / 7\ / \ v{2n/d) (pinid) , . , 

c{j,2n/d) = f^{t2n/d,j)—r. 7 = -fJ-[tn/d,j)—r. r = -c{j,n/d). 

"Pihn/dj) fKtn/dj) 

Now suppose that j G 2N. We have gcd{2n/d,j) = 2 gcd{n/d, j) and also t2n/d,j = ^n/dj- This yields 
directly that c{j,2n/d) = c{j,n/d). □ 



Theorem 11 For an arbitrary integral circulant graph ICGn{D) there exists an odd number < j < n — 1 
with Xj also odd if and only if there is a divisor d & D satisfying d/2 ^ D and 2d ^ D. 

Proof. 

(=^>:) Suppose that for every ci G D it holds that either 2d e D ov d/2 G D. It follows that D = DiU2Di. Let 
J G 2N + 1. According to Lemma [TOl we have c{j, n/d) = —c{j, 2n/d) for any d G D such that n/d G 2N + 1 
and therefore 

= Yl ^/^) = Yl "/^) + 2"/^) = G 2N. 

deDiU2Di deDi 
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We conclude that all eigenvalues with odd indices are even. 

(■^0 Let D'l = {deD\n/d£m + 2=f- 2d £ D} and D" = D'l U 21)'/. 

Let D' = D \ D". By the assumption, D' is a non-empty set. According to Lemma [10] it holds that 
c{j, n/d) + c{j, 2n/d) G 2N for every d G D" such that n/d £ 2N + 1. Denote by d^a^. = maxD'. Since d^aa; 
is a divisor of n, it can be represented in the form dj^^x = 2^"^^^ • • • p^*" where < /3j < for i = 1, . . . ,k. 
Without loss of generality, we can suppose that there exists 1 < s < A: such that Pi < Ui for i = 1, . . . , s and 
/3j = Oj for z = s + 1, . . . , /c. Then we can write n/d'^^^ = 2°'°~^^p1^~^^ ■ ■ ■ ps"~^'' ■ Denote by 

ai — Bi — l a,— 0,-1 Os+i oi. 

It holds trivially that < jo < n — 1. Lemma [9] yields directly that c{jo,n/d'j^a^) is odd since n/d'^^^ > 2. 

Suppose that D'\{d'j^^^} ^ and let d E D'\{d'^^^} be an arbitrary divisor with its prime factorization 
d = 2'yopf ...pf. 

We will show that there exists 1 < i < k such that < 7^ < /3j < a^. Suppose this is not the case, which 
means that < /3j < 7^ < aj for 1 < i < k. If /3o < 70 then d^a^ | d, which is a contradiction. Similarly, 
if /3o = 7o + 1 then it holds that (imaa; I 2c?, which implies d = d'^g^^/2, providing a contradiction with the 
definition of the set D' . 

Let i be an arbitrary index such that 7, < /3j. Suppose that Oi — /3i > 1. Then i < s and Sp.{n/d) = 
Oii — 2. Since Sp^{jo) = — /3j — 1 > Oj — 7^ — 1 = Sp^{n/d) — 1, from Lemma Owe can conclude that 
c{jo,n/d) is even. 

Now suppose that Oi = /3j. Then i > s and Sp.{n/d) = — 7i > 1. Again since Sp.{jo) = Oj > 
aj — 7i — 1 = Sp^{n/d) — 1, Lemma [9] yields that c{jo,n/d) is even. 

This implies that there is an odd index Jq such that c{jo,n/d'^g^^) is odd and c{jo,n/d) is even for every 
d £ D' \ {d'j^a^}. Now we have 

Ajo = c{jo,n/d'^^.J + ^ c{jo,n/d) + ^ c(jo, n/d) G 2N + 1, 

since both sums in the last expression are even. 

If dmax is the only divisor contained in D, the above sum is reduced to 

Ajo = c{h,n/d'^^^) + c{jo,n/d) G 2N + 1, 

deD" 

and it is still even. □ 



Let us mention two direct consequences of the previous theorem. The first one is actually the contra- 
positive of the assertion of the theorem. 

Lemma 12 All eigenvalues o/ICG„(D) on odd positions are even if and only if D = Di U 2Di. 

Lemma 13 Let n/2 G D. All eigenvalues of lCGn{D) on odd positions are odd if and only if D = D* U 
2DI U {n/2} where Dl = Di\ {n/2}. 

Proof. Suppose that all eigenvalues of ICG„(-D) on odd positions are odd. Denote by Xj the eigenvalues 
of the graph. Now, consider the integral circulant graph ICG„(-D') where D' = D \ {n/2}. Denote by 
the eigenvalues of integral circulant graph ICG„(L''). Since 

'''^-^^)-\ 1, 2\j ' 1, 2\j ■ (1^) 



it holds that 
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This gives that all eigenvalues of ICG(Z)') on odd positions are even. According to Lemma [T2] we have 
that D' = DlU 2DI where Dl = Di\{n/2} which completes the first part of the proof. 

The converse of the assertion can be proven analogously. Let D = Djj'U2L'|U{n/2} and D' = D\{n/2}. 
According to Lemma [12] all eigenvalues of ICG„(L'') on odd positions are even. Furthermore, using the 
relation ()14p we have that all eigenvalues of ICG„(L') on odd positions are odd. □ 

4.2 Perfect quantum distance, spectrum and characterization of ICG„(D) having PST 

The main result of this section is characterization of ICGn(-D) having PST. In addition, we also calculate 
perfect quantum distance of integral circulants having PST. Thus, according to Theorem [8] from now on we 
assume that n is divisible by four. 

Theorem 14 Let ICGn(D) have PST. Ifn/2 G D then all eigenvalues are odd, otherwise they are even. 
Moreover, in the first case the eigenvalues on odd positions are equal to —1 and in the second one the 
eigenvalues on odd positions are equal to 0. 

Proof. According to Proposition [2] we have Ai = X^deD ^(^/^)' Since 4 | n/d for d £ D2 U D-^ we 
conclude that iJ,{n/d) = and therefore Ai = J2deDi Using Proposition [2] once again we see that 

^2 = EdeDo Kn/d) + T^deDi Kn/^d) + EdgDjuSg 2^(^/2^). For d e D3 we have 4 | n/2d, which yields 

A2 - Ai = {ii{n/2d) - nin/d)) + 2 ^ n{n/2d) G 2N. 

By Lemma El all the differences A^+i - Ai G 2N for < z < n - 2, since lCGn{D) has PST. If n/2 ^ D 
then Ao G 2N and thus all the eigenvalues are even, otherwise Aq G 2A^ + 1 and all the eigenvalues are odd. 

Let j G 2N + 1. For d £ D2U D3 we have that 4 | tri/dj thus c(j, n/d) =0. This reduces the formula 
for the j-th eigenvalue to 

Xj = J2 c(.j,n/d). 
der)i 

The condition n/2 ^ D yields that Aq G 2N and according to the first part of the proof all the eigenvalues 
are even. 

Let //j, < j < n — 1, be the eigenvalues of the integral circulant graph ICG„(-Di). Thus Aj = for 
any odd < j < n — 1, so /ij G 2N for j G 2N + 1. From Lemma [12] we conclude that all the eigenvalues 

on odd positions are even if and only if Dq = 2Di and Di = DiU 2Di. By the first part of the proof of 
Theorem [11] we deduce that fij = for j G 2N + 1 and consequently Xj = for j G 2N + 1 . 

Analogously, if n/2 G D we obtain that Xj G 2N + 1 for < j < n — 1 and so G 2N + 1 for odd 
< j < n - 1. Now, Lemma [13] yields D = D^U 2DI U {n/2} where D*^ = Di\ {n/2}. 

Now, consider the integral circulant graph ICG„(D') where D' = D \ {n/2} = D\ U 21)^. Denote by 
A^- the eigenvalues of the integral circulant graph ICG„(L''). From we obtain that Xj = X'^ — 1 for 
J G 2N + 1. But, according to the first part of the proof we have that A^ = for j G 2N + 1 and so Aj = — 1 
for j G 2N + 1. 

□ 

From the proof of the last theorem we can derive the following important corollary 

Corollary 15 //ICG„(L») has PST then Dq = 2{Di \ {n/2}). 

Using Theorem [T3] we can establish a more precise criterion for the characterization of integral circulant 
graphs having PST, than the one given by Corollary [7] 

Lemma 16 ICGn(-D) has PST if and only if there exists an integer k > 1 such that one of the following 
conditions are satisfied 

i) 5*2 (A2j) = k and X2j+i = 0, if n/2^ D 
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ii) S2{X2j + I) = k and X2j+i = —1, if n/2 S D 
for 0<j< n/2. 

Lemma 17 Let n be an even number, d be a divisor of n and ni = n/2. For an even number < j < n — 1 
the following equalities hold 

1. c{j,n/d) = c{3/2,^), ifdGDo 

2. c{j,n/d) = c{j/2,ni/d), if d e Di 

3. c{j,n/d) = 2c{j/2,ni/d), ifdeD2UD3. 
Proof. 

1. Suppose that d £ Dq. Then gcd(n/(i,j) = gcd(^,j72) and 

n 2ni 
" dgcd(nMi) " 2f gcd(J^,j72) " 

Furthermore, it holds that ip{n/d) = V'(^) and so c{j,n/d) = c{j /2, 

2. Suppose now d £ Di. Then gcd(n/(i, j) = 2gcd{ni/d,j/2) and 

n 2ni 
" (igcd(n/d,i) " d2gcd(ni/(i,i/2) " 

We also conclude that (p{n/d) = (p{2^) = ip{ni/d) and so c{n/d,j) = c{ni/d, j /2). 

3. Suppose now d £ D2U D3. By the preceding case we obtain that gcd{n/d, j) = 2gcd{ni/d,j/2) and 
tn/d,j = ^ni/d,j/2- Let k = S2{n/d) and n' be an integer such n/d = 2^ n' /d. Notice that A; > 2, since 
d G D2UD3. Also, n'/d is odd and thus gcd{2'' , n' / d) = 1 which yields ifin/d) = ip{2^''^) = 2'^-V(^)- 
Furthermore we have 

^{n/d) = 2'=-V(^) = 2^{2^-^)^C^) = 2^(2'^-'^) = 2^{m/d) 
and so c{j,n/d) = 2c{j /2,ni/d), ii d G D2U D3. 

□ 



Lemma 18 //ICG„(L») has PST then Di = 2{D2 \ {n/4}). 
Proof. 

Let \j be an eigenvalue of ICGn(-D) where < j < n — 1. 
According to Corollary [T^] it holds that Dq = 2{Di \ {n/2}) and so 

Xj = ^ c{j,n/d) + ^ c{j,n/d) + ^ c{j,n/d). 
de2{Di\{n/2}) rfeDi dGDaUAi 

From the relation ()13p we have c(j, 2) = 1 for j G 2N. 

For d e 2{Di \ {n/2}), let d = 2d' where d' e Di \ {n/2}. Then c{j,n/d) = c{j,ni/d') and according to 
Lemma [T71 f part 1.), we have c{j,ni/d') = c{j /2, ^77^) = c{j /2,nl/d'). 

II d £ Di using Lemma [T71 (part 2.) it holds that c{j, n/d) = c(j/2, nl/d). 

Finally, if d £ D2LS using Lemma [T71 (part 3.) it holds that c{j,n/d) = 2c{j /2,nl/d). 

Taking the discussion above into account we obtain 
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^ r cU2)+2E,,,,,\{„/2}uD.uD3^(^V2,ni/d) = 2A;./2 + l, n/2eD 
' 1 2E,eDAW2}uD.u53<jV2,ni/d)=2A;./„ n/2^D ^ 

where ni = n/2 and A^ are the eigenvalues of the integral circulant graph ICGnj(L'') where D' = {Di \ 
{n/2}) U L>2 U I^s. _ 

We conclude that gcd(n/(i, j) = 2 for j G 4N + 2 and d G Ds, hence that 4 | tn/d,ji aiid finally that 
c{j,n/d) = 0. This yields that 

>^ >^ f 2A'.,„ + 1, n/2 e D 

where A^- = J2deDi\{n/2} ^{ji'^i/d) + l^cieD2 '^i/d). This means that the eigenvalues A^- for odd < j < 
ni — 1 coincide with the eigenvalues on odd positions of ICG„j(Z)i U \ {'^/2})- Let us denote by /ij the 
eigenvalues of ICGn^Di U L>2 \ {n/2}), for < j < ni - 1. 

Prom Lemma [16] it follows that S2{\j) = k 01 S2{Xj + 1) = k, for j G 4N + 2 and some integer A; > 1, 
depending on whether n/2 ^ D or n/2 G D. Furthermore, we have that either 52(A^y2) = — 1 or 
S2{^'jl2 + ^) = k — 1 ior odd < j /2 < ni — 1. Since = X'j for odd < j < ni — 1 (eigenvalues /Xj on odd 
positions have the same parity), Lemma [T2] and Lemma [13] yield that Di = 2{D2 \{ni/2}), which completes 
the proof. 

□ 



Theorem 19 //ICG„(L>) has PST then either n/2 e D or n/4 G D. 
Proof. 

Case 1. n/2 D. Suppose also that n/4 D. Since n/2 ^ D, according to Corollarv 1151 it holds that 
Dq = 2Di and using Proposition [2] we have A2 = Ede2Di /^(^/^^) + EdeDi l^{n/2d) + 2J2d(zD2UDs n{n/2d). 
For d £ D3 we conclude that 4 | n/2d and so X^^eDa f^i^/'^^) ~ 0- Now, the formula for the eigenvalue A2 
becomes A2 = "^deDi ^(n/2(i) + J2deDi /^(n/2(i) + 2J2deD2 /^(^/2d)- Since n/A ^ D, Lemma [TH] now leads 
to 

A2 = 2( ^ ^x{n/2d) + ^in/2d)) = 2( ^ ^(n/4(i) + ^ ^(n/2d)) = 0. 

de2D2 d&D2 d&D2 d&D2 

Case 2. n/2 G Z?. Suppose also that n/4 G -D. This gives Z^o = 2{Di \ {n/2]) and Di = 2{D2 \ {n/4}), 
which follows from Corollarv 1151 and Lemma [18] In this case A2 can be written as follows 

X2=Yfi{n/d)+ fi{n/2d) + 2 ^ /i(n/2d) + ^(1) + 2/i(2) = + 1 - 2 = -1. 

deD[ d(^Di\{n/2} deD2\{n/4:} 

Notice that EdeA) ^^('^/^) + EdeDi\{n/2} f^{n/2d) + 2 EdGD2\{n/4} = according to Case 1 of the proof. 
In both cases we conclude Ai = A2 (according to Lemma [T6]l and this is a contradiction according to 

Corollary [5] □ 



Perfect quantum communication distance (PQCD) of an arbitrary pair of vertices a and b is the distance 
d{a, b) if a perfect state transfer exists between them. If we consider a circulant network with identical 
couplings PST occurs only between vertices b and 6 + n/2 for 0<6<n/2 — 1 (Theorem H]) . For the integral 
circulant graph ICG„(D) PQCD of b and b + n/2 is equal to one, if n/2 G D. Otherwise, we have that 
n/4 £ D (Theorem I19p and thus the path b,b + n/4:,b + n/2 shows that PQCD is equal to two. In both 
cases PQCD is independent of the order of the graph. 

Now, we are ready to describe the spectrum of integral circulant graphs. The criterion for existence of 
PST in integral circulant graphs that we will use in the next two theorems is given by the following lemma. 
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Lemma 20 ICG„(D) has PST if and only if one of the following conditions holds 
i) \2j G 4N + 2 and A2j+i = 0, ifn/2 D 
ii) \2j G 4N + 1 and \2j+i = -1, ifn/2 G D 
for 0<j< n/2. 

Proof. In the proof we will use the same notation as in the proof of Lemma [THl Suppose that ICGn(D) 
has PST. 

From Lemma [16] it follows that \j = or Aj = — 1 depending on whether n/2 ^ D or n/2 G D, for any 
odd < i < n - 1. 

Now suppose that n/4 G D. According to Theorem [19] we have that n/2 ^ D. This implies that Aq G 2N 
and hence all the eigenvalues of ICG„(L') are even for j G 2N, which follows from Lemma [TU] Now, we 
proceed using the proof of Lemma [18] We have that 

= 2Aj/2 

where A^ are the eigenvalues of the integral circulant graph ICGn^(D') where D' = Di U D2 U D3 \ {n/2} 
and ni = n/2. 

Since n/4 = ni/2 e D it follows that Aq G 2N + 1 which further implies that Aq G 4N + 2. Finally, from 
Lemma [16] we conclude that Xj G 4N + 2 for j G 2N. 

If n/2 £ D then Aq G 2N + 1 and hence Xj G 2N + 1 for j G 2N. This further implies that Xj = 2X'j^^ + 1. 
According to Theorem 1191 it holds that n/4 = ni/2 D and so Ag G 2N. This implies that Aq G 4N + 1. 
Finally, from Lemma [16] we conclude that Xj G 4N + 1 for j G 2N. 

If any of i) or ii) holds, it can easily be seen that Aj+i — Aj G 4N + 2 which implies that ICG„(D) has 
PST, according to Lemma [U] 

□ 



Theorem 21 Let D be a set of divisors of n such that n/2, n/4 D. Then ICG„(Z) U {n/4}) has PST if 
and only i/ICG„(Z) U {n/2}) has PST. 

Proof. 

Let Xj, jdj and be the eigenvalues of graphs ICG„(D U {n/4}), ICG„(D U {n/2}) and ICG„(D), 
respectively. We have the following relations between these eigenvalues: Xj = Uj+c^j, 4) and Hj = Vj+c{j, 2). 
This yields that ^ij = Xj - c{j, 4) + c{j, 2) for < j < n - 1. 



By direct computation we show that 

4 



td. 



gcd(4,j) 
From this it follows that 




Hj 



c(j,4) 



' Aj -l, JG2N + 1 
Xj + 3, j G 4N + 2 
A,-l, jG4N 



0, j G 2N + 1 
-2, j G 4N + 2 . 
2, j G 4N 



(17) 



(18) 



The following two facts can now be easily deduced: for j G 2N +1, Xj = if and only if fij = —1 and for 
j G 2N, Aj G 4N + 2 if and only if fij G 4N + 1. To complete the proof it only remains to apply Lemma (201 

□ 



Finally we can state the next of our main results. 

Theorem 22 ICG„(Z)) has PST if and only if n e 4N, = 2D^, Dq 
n/2 G D, where = L»2 \ {n/4} and = Di\ {n/2}. 



AD 2 and either n/A G D or 
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Proof. 

(=^:) This is an easy consequence of Theorem [HI Lemma [TSl Coroharv 1151 and Theorem I19[ 

(<^:) According to Theorem 12 H this impUcation is sufficient to prove for n/4 G D. Furthermore, by Theorem 

[I9]we have n/2 D. 

Let < J < n — 1 be an odd number. For d £ D2 U Ds, we conclude that c{j,n/d) = 0, which fohows 
from the fact that 4 | tn/d,j- This imphes that 

= Yl c{j,n/d) + c{j,n/d) = c{j,n/2d) + c{j,n/ d) = 0. 

The last equality follows from Lemma [TOl 

Let < j < n — 1 be an even number. We have 

^3 = 4j^n/d) + ^ c{j,n/d) + c(j, n/d) + c(j, 4) + ^ c{j,n/d). 

de2Di deDi d£D2\{n/A} deDa 

From Lemma [10] and relation (|17p it follows that 

A,- = 2 ^ c{j, n/d) + '^Z^) + Yl n/d) ± 2. 

Now using Lemma [T7| we obtain 

A,- = 2 ^ c(i/2,ni/(i) + 2 ^ c(j72,ni/(i) + 2 ^ c(j72,ni/d) ± 2 

dGDi=2(D2\{n/4}) deD2\{n/4} dggg 

= 2 (c(j/2, + c(j72, ni/2d)) + 2 c(j/2, m/d) ± 2, 

rfe-D2\{n/4} 

where nl = n/2. 

Let j e 4N + 2. By Lemma [JO] we obtain that c(j/2, ni/d) + c(j72, ni/2(i) = 0, since j/2 e 2N + 1. For 
d £ we conclude that 4 | tni/d,j72 aiid c{j /2,ni/d) = 0. Finally, we conclude that Xj = c(j, 4) = —2. 
If j E 4N, according to Lemma [10] we show that 

Ai = 4 c(i/2,ni/<i) + 2 c(j72,ni/d) + 2. 

d€D2\{n/4} deD3 

Furthermore, using Lemma [T7| we have c{j /2,ni/d) = 2c{j /4:,ni/2d). In either case we conclude that 
\j S 4N + 2 for j G 2N. Now, direct application of Lemma [20] completes the proof. 

□ 

According to the previous result, we notice that graph ICG„(-Dn \ {n/2}) has PST for n G 4N. This class 
of circulant graphs is known as cocktail-party graphs (see |16]). 

5 Conclusion 

In this paper we continue to address the question of when circulant graphs can have perfect state transfer, and 
improve the necessary condition of n being divisible by 4 given in [17] , Theorem [22] completely characterizes 
the graphs ICG„(D) having PST. This result includes the classes of graphs having PST found in [3] [T5] [T7]. 

From the above characterization we can calculate the number of integral circulant graphs of a given 
order having PST. If n E 8N, by the rule of product, the number is equal to the product of the cardinalities 
of the power sets of {d : d | n, n/d G 8N} and {d : d \ n, n/d G 8N + 4} \ {n/4} times two. If n G 8N + 4 
the number is equal to the cardinality of the power set oi {d : d \ n, n/d G 8N + 4} \ {n/4} times two. In 
either of cases we have two possibilities since either n/2 G D or n/4 G D. Thus, for a given number n the 
number of integral circulant graphs ICGn(-D) having PST is given by the following formula 
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2^(t), nG8N + 4 
n G 



|ICG„(Z?)| = ^.(tM^); _ _ , (19) 



where T(n) denotes the number of the divisors of n. 

We can see from the formula that for some values of n E 8N (for example n = 96, 120, 144, 160, 168, 192, . . .) 
there is a great number of graphs having PST, while for some other values n G 8N + 4, there are only 2 such 
graphs. However, the number of ICG„(-D) having PST is asymptotically equal to the number of ICG„(D) 
of a given order n. The last conclusion follows from Corollary 7.2 given in [20], where it was shown that 
there are at most 2"^^"^"^ integral circulant graphs on n vertices. 

It is worth mentioning that the maximum value of the perfect quantum communication distance (i.e. the 
maximal distance between vertices where a perfect state transfer occurs) is equal to 2 for every n S 4N. 
Thus, it is still an open problem whether one can construct a network with identical couplings in which 
any quantum state can be perfectly transferred over a larger distance than 21og3n, obtained in [3 [8] for 
two-link hypercubes with n vertices. 

An improvement of the perfect quantum communication distance is made in [7j by considering fixed 
but different nearest-neighbor couplings. A similar approach used on circulant graphs (having a weighted 
adjacency matrix) might also enlarge the perfect quantum communication distance. Many recent papers 
propose such an approach |1 H 114 1 [151 116] . First results concerning characterization and finding new classes of 
weighted circulant graphs are given in [5] . Characterization of integral circulant graphs (moreover circulant 
graphs) having PST is the first step in describing the more general class of weighted integral circulant graphs 
having PST. The approach to this problem should, in our opinion, use the interplay of graph and number 
theory. 
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